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ROBOT  MANIPULATOR  CONTROL 

Claude  Samson  and  J.M.  Ibarra  Zannatha,  Institut  de  Recherche  kn 
Informatique  et  Systemes  Aleatoires,  Rennes,  France,  have  developed  a  gen¬ 
eral  approach  to  the  problem  of  controlling  robot  manipulators.  Samson  and 
Zannatha  have  separated  the  larger  problem  into  two  issues:  controlling  a 
known  time-invariant  linear  system  and  modelling. 

The  two  sub-problems  are  related  when  control  robustness  is  considered 
with  respect  to  modelling  errors.  Samson  and  Zannatha's  approach  allows 
derivation  of  several  adaptive  or  nonadaptive  control  schemes  described  in  the 
literature  and  suggests  new  schemes.  It  applies  to  continuous  and  discrete 
cases  and  shows  possibilities  for  computing  controls  from  a  set  of  operational 
coordinates  of  the  manipulator. 

This  report  is  a  translated  and  abbreviated  version  of  a  paper  titled 
"Sur  La  Commande  Des  Manipulator"  by  Samson  and  Zannatha. 


Mathematical  Model  of  a  Manipulator  With  N  Degrees  of  Freedom 

We  begin  with  the  general  form  of  a  mathematical  model  of  a  manipulator: 

M(x)  x  +  V(x,x)  =  u  (2.1) 

in  which  x  is  the  vector  of  n  generalized  coordinates,  x  is  its  time  derivative 

(speed  vector),  x  is  its  second  derivative  (acceleration  vector),  and  M(x)  is 
the  kinetic  energy  or  inertia  matrix,  which  is  always  positive  definite. 
V(x,x)  is  the  vector  combining  the  Coriolis  and  centrifugal  forces;  u  is  the 
force  and  moments  vector  acting  on  the  manipulator,  other  than  V(x,x).  One 
can  decompose  u  into  the  sum  of  three  terms: 


u  =  u  (x)  +  u,(x)  +  u 
g  t  a 


(2.2) 


where  u  (x)  is  the  force  of  gravity,  uf(x)  is  friction,  and  u  represents 
9  T  a 

forces  produced  by  the  motors  on  the  movement  of  the  robot.  The  force  of 

interest  here  is  u  . 

a 

If  we  combine  (2.1)  and  (2.2),  the  manipulator  model  can  be  written: 


M(x)  x+N (x,x)  =  u. 


(2.3) 


where 


N(x,x)  *  V(x,x)  -  u  (x)  -  u.(x).  (2.4) 

9  f 

Dividing  the  Initial  Control  Problem  Into  Linear  Control  and  Modelling 
Problems 

Equation  2.3  is  replaced  by  the  following  system  of  equations  equivalent 

to  it. 
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If  M(x)  and  N(x,x)  are  known  or  easily  calculated  at  each  instant,  then 
it  is  sufficient  to  calculate  the  control  u'  associated  with  an  invariant  linear 
system  of  second  order  (2.5a),  and  then  to  calculate  u  using  (2.5b,  c,  and 
d).  a 

If  M(x)  and  N(x,k)  are  not  known  and  are  too  complex  to  calculate  at 

each  instant,  we  can  work  with  approximate  matrices  M  and  N  and  use  the 
following  working  model  derived  from  the  mathematical  mode)  (2.5): 


u  =  M  u’ 
al 

Ua2  =  N 

Ua  =  Ual  +  Ua2* 


(2.6a) 
(2,6b) 
(2.6c) 
(2. 6d) 


Thus  the  problem  of  control  of  the  manipulator  is  the  sum  of  two  sub- 
problems: 

1.  A  linear  control  problem  associated  with  equation  (2.6a),  in  which  we 
can  use  the  classic  theory  of  linear  and  invariant  systems  control  to  compute 

u'. 

2.  A  modelling  problem  associated  with  computing  the  matrices  M  and  N. 

The  control  will  be  only  as  accurate  as  the  approximation  of  the  matrices  M 
A  • 
and  N  to  the  matrices  M(x)  and  N(x,x). 

Although  the  two  problems  require  different  techniques,  they  are  not 
totally  independent.  Errors  associated  with  the  modelling  problem  may  have 
more  or  less  serious  consequences  depending  on  the  type  of  control  chosen 
for  u'. 

Two  Strategies  of  Control  for  Following  a  Reference  Trajectory:  Consigned 


Control  and  Control  With  Reference  Model 

The  strategy  for  control  u  is  determined  by  the  computation  of  u' .  To 

3 

simplify  the  study  we  recall  the  principle  of  two  methods  most  often  used  for 
the  control  of  manipulators:  consigned  control  and  control  with  reference 
model. 

Consigned  Control 

When  the  manipulator  is  first  at  rest  [x(t  )  =  o] ,  we  wish  to  have  x 

0 


return  to  a  consigned  x  according  to  a  certain  dynamic.  We  consider  the 
following  control  u': 

u'  =  -Kp(x-xc)  -  Kvx  (4.1) 

(Kp  and  Kv  are  diagonal  matrices). 


We  obtain,  following  (2.6a),  the  transfer  function: 

x(s)  =  (K  /  (s2+K  s+K  ))x  (s).  (4.2) 

p  v  p  c 

The  time  constant  and  the  overshoot  of  the  indexed  response  of  the 
system  (4.2)  are  regulated  by  the  choice  of  gains  Kp  and  K^.  The  gains  are 

therefore  entirely  determined  by  the  dynamic  of  the  return  to  the  consigned 
value. 

Control  With  Reference  Model 

The  method  consists  of  giving  the  equation  of  a  reference  system  (or 
reference  model)  of  the  same  dimensions  as  the  model  of  the  manipulator,  and 
computing  the  control  u*  so  as  to  cancel  the  error  e  between  the  position 
vector  x  of  the  manipulator  and  the  position  vector  x^  of  the  reference  model. 

The  choice  of  the  reference  model  must  be  compatible  with  the  structure 
of  the  manipulator.  There  should  be  a  bounded  control  u'  allowing  cancel¬ 
lation  of  the  error  e  (see  H.  Erzberger,  "Analysis  and  Design  of  Model  Fol¬ 
lowing  Control  Systems  by  State-Space  Techniques"  Proceedings  of  the  Joint 
Automatic  Control  Conference  (1968],  pp  572-581).  Referring  to  equation 
(2.6a),  an  evident  choice  for  the  reference  model  is: 

x  =  u  .  (4.3) 

r  r 

Let 

e  =  x  -  x  (4.4) 

r 

subtracting  (4.3)  from  (2.6a)  we  obtain 

e  =  u *  -  ur.  (4.5) 

A  possible  choice  for  the  control  u'  is  then: 

u '  =  u  -  Le-Le  (4.6) 

r  p  v 

where  L  and  Ly  are  diagonal  matrices.  The  dynamic  equation  of  the  error 
e  is  givfen  by: 
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(4.7) 


Using  the  Routb-Hurwitz  criteria,  it  can  be  shown  that  (4.7)  is  the  equation 
of  an  exponentially  stable  system,  and  the  elements  of  the  diagonal  matrices 
L  and  are  positive.  The  choice  of  higher  gains  leads  to  faster  exponen¬ 
tial  convergence  toward  zero. 

The  control  in  the  reference  model  acts  uniquely  as  a  regulator  because 
its  role  is  to  bring  the  error  to  zero.  The  reference  control  ur  is  required 

to  follow  a  reference  trajectory.  Unlike  the  consigned  control  method,  the 
choice  of  gains  and  are  subject  to  no  constraints  coming  from  the 

dynamic  of  the  reference  trajectory.  It  can  be  shown  that  the  method  of 
consigned  control  is  a  special  case  of  the  reference  model  method  obtained  by 
choosing 


u 


r 


with 


K  x 


P  r 


K  x 


v  r 


+  K 

p 


x 

c 


Control  Schemes  as  a  Function  of  Modelling 

The  relations  (2.6b,  c,  ai.d  d)  and  (4.1)  can  be  combined  to  obtain  the 
consigned  control: 

u  =  -MK  (x-x  )  -  MK  x  +  N.  (5.1) 

cl  p  C  V 

The  relation  (2.6b,  c,  and  d)  with  (4.6)  give  the  reference  model  con¬ 
trol: 


u  *  -  ML  e  -  ML  c+  Mu  +  N.  (5.2) 

a  p  v  r 


First  Modelling 

The  classic  controls  with  constant  coefficients  are  those  computed  from 
simplified  models  of  manipulators  using  the  following  hypothesis: 

A 

Hj  :  M(x)  =  M  :  constant  diagonal  matrix 


H 


2 


H 


3 


ug(x)  a  0;  uf(x) 
V(x,x)  s  0 


1 


N 
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Hypothesis  is  justified  for  small  displacements  and  when  the  N  de¬ 
grees  of  freedom  of  the  manipulator  are  practically  decoupled.  It  also  sup¬ 
poses  that  the  matrix  of  inertia  is  only  slightly  sensitive  to  the  variations  of 
mass  carried.  Hypothesis  indicates  that  the  gravity  terms  (often  compen¬ 
sated  in  part  by  counterweights),  and  the  friction  terms  ^are  negligible. 
Finally,  hypothesis  H^  is  justified  if  the  displacement  speed  x  of  the  manip¬ 
ulator  is  small. 

Associating  this  model  with  (5.2)  gives  a  nonclassical  scheme  of  model 
reference  control  with  constant  coefficients. 

Ceocnd  Modelling 

Using  the  mathematical  model  of  the  manipulator  we  calculate  for  each 
point  of  the  trajectory  the  matrices  M(x)  and  N(x,x),  which  gives: 

M  =  M  (x) 

N  =  N (x,x)  . 


Controls  using  the  model  are  called  dynamic  controls;  it  is  possible  to 
have  dynamic  control^  with  a  reference  model.  The  computation  of  the  mat¬ 
rices  M(x)  and  N(x,x)  is  complex  and  in  some  situations  cannot  be  done — for 
example,  when  the  mass  and  shape  of  the  object  carried  vary  greatly  and  are 
unknown. 

Third  Modelling 

We  identify  on-line  the  matrices  M  and  N  beginning  with  a  parameter¬ 
ization  of  the  matrices  established  in  advance.  The  controls  thus  synthesized 
are  adaptive.  More  generally,  we  can  design  by  adaptive  control  all  controls 
requiring  measurements  made  on  the  manipulator  to  identify  its  structure  or 
the  control  law  (see,  for  example,  S.  Dubowsky  and  D.J.  Des  Forges,  "The 
Application  of  Model -Referenced  Adaptive  Control-  to  Robotic  Manipulators," 
Journal  of  Dynamic  Systems,  Measurement  and  Control,  Vol  101  [September 
1979],  pp  193-200).  According  to  the  definition,  the  dynamic  controls  previ¬ 
ously  discussed  are  not  adaptive.  The  mathematical  model  of  the  manipulator 
from  which  they  are  computed  will  be  a  function  of  the  measurements  of  x  and 
x. 

The  idea  of  adaptive  control  is  in  general  motivated  by  two  considera¬ 
tions: 

1 .  No  mathematical  model  of  the  manipulator  completely  reflects  the 
reality  because  of  the  difficulty  of  modelling  terms  (such  as  friction)  and 
changes  in  the  course  of  use  (for  example,  picking  up  a  load). 

2.  An  exact  model  of  the  manipulator  is  less  important  for  control  than 
an  accurate  identification  of  the  structure  of  the  manipulator  arising  from  its 
principal  characteristics. 

There  is  no  general  method  for  on-line  identification  of  M  and  N.  There 
are,  however,  techniques  for  estimating  recursively  a  vector  e  of  constant 
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parameters,  or  slowly  changing  variables  when  we  have  a  relation  of  the 
following  form  (see  L.  Ljung,  and  T.  Soderstrom,  Theory  and  Practice  of 
Recursive  Identification  (Cambridge,  MA:  MIT  Press,  1 982 1 ) : 

y  (t)  -  <J>  ( t )  T0  ,  (5.3) 

where  y(t)  is  a  known  vector,  and  4>(t)^  is  a  known  matrix  of  observation 
with  dimensions  compatible  with  those  of  y ( t )  and  6. 

Such  a  relation  will  normally  be  obtained  after  parameterization  of  the 
dynamic  equation  of  the  manipulator: 

M  x  +  N  =  u  ,  (5.4) 

cL 

which  leads  to  the  adaptive  control  schemes  called  "indirect."  We  now  give 
several  examples  of  parameterization. 

Example  1 .  The  vector  0  is  made  up  of  the  totality  of  components  of  the 
matrices  M  and  N.  Relation  (5.4)  can  be  written  in  the  form  (5.3)  with  y(t) 

=  u  (t)  and  $  (t)^  made  up  of  unity  and  the  components  of  acceleration  x(t). 

3 

In  the  continuous  case,  the  parameterization  assumes  knowledge  of  the  accel¬ 
eration  x(t).  In  using  the  symmetry  of  the  matrix  M,  the  number  of  compon¬ 
ents  entering  into  0  can  be  reduced  substantially.  Also,  the  knowledge  of 
zero  components  of  M,  from  consideration  of  the  geometry  of  the  manipulator, 
reduces  the  size  of  0. 

With  such  parameterization,  the  matrices  M ( t )  and  N(t)  are  directly 

A 

obtained  from  the  identified  vector  0(t). 


Example  2.  If  the  gravity  and  friction  forces  are  negligible,  use  of  (2.4) 

reduces  N(x,x)  to  V(x,x),  and  we  can  write  the  following  (see  R.  Horowitz 
and  M.  Tomizuka,  "An  Adaptive  Control  Scheme  for  Mechanical  Manipulators 
Compensation  of  Nonlinearity  and  Decoupling  Control,"  American  Society  of 
Mechanical  Engineers  (ASME),  Paper  No.  80,  Wa/DSC-6  (1980]): 


(5.5) 


where  the  matrices  v'(x)  are  symmetric  and  only  depend  on  the  position  x. 

A  possible  parameterization  is  then  obtained  by  putting  into  the  vector  0 

the  unknown  components  of  M  and  of  the  matrices  v'(x).  The  parameu  r- 
ization  leads  to  a  vector  0,  generally  larger  than  in  the  preceding  example. 
Its  advantage  is  that  the  components  entering  in  0  are  not  functions  of  the 
speed  x,  and  for  movements  of  small  amplitude  we  can  expect  the  vector 
6(x)  to  be  nearly  constant. 
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Example  3.  For  a  given  load,  the  parameters  in  the  mathematical  model  of 
the  manipulator  are  fixed  and  we  can  try  to  identify  them  on-line  and  enter 
them  into  the  vector  6.  For  a  given  configuration,  e  is  constant  and  thus 
independent  of  x  and  x  and  the  computed  control. 

Identifying  6  in  the  above  case  involves  computation  at  each  instant  of 
the  coefficients  multiplying  the  components  of  0  in  the  mathematical  model  of 
the  manipulator.  The  computed  controls  are  dynamic  adaptive  controls,  which 
in  some  cases  are  too  complex  for  easy  computation. 

Observation.  The  set  of  adaptive  controls  that  we  have  considered  are 
studied  in  more  detail  in  J.M.  Ibarra  Zannatha,  "Sur  La  Command  Des  Robtos 
Manipulateurs"  (Doctoral  thesis,  Univ.  of  Rennes,  April  1982).  The  following 
observations  are  based  on  simulations  during  which  many  identification 
algorithms  were  tested. 

1 .  The  parameterizations  in  examples  1  and  2  do  not  permit  easy  identi¬ 
fication  of  M  and  N  (poorly  or  noisily  followed).  The  poor  identification  is 
due  to  two  main  causes:  the  parameters  that  we  seek  to  identify  vary  too 
rapidly  with  rapid  movement  of  the  manipulator,  and  the  parameters  depend 
on  the  intermediate  values  of  x  and  x.  With  the  adaptive  control  depending 
on  estimated  parameters,  the  identification  problem  for  the  two  parameter¬ 
izations  is  ill  conditioned. 

2.  In  spite  of  the  mediocre  performance  of  the  identification,  the  adap¬ 
tive  controls  with  the  model  of  reference  retain  good  performance  when  the 
gains  and  are  large  (the  reference  trajectory  is  closely  followed).  The 

interpretation  of  the  phenomenon  (observable  in  many  existing  studies)  is 
important  to  understanding  the  problems  in  the  control  of  manipulators.  We 
attribute  this  largely  to  the  robust  character  of  the  large  gain  control. 

Robustness  of  Large  Cain  Control 

A  rapid  heuristic  study  justifies  the  robust  character  of  large  gain 
controls  with  respect  to  errors  in  the  modelling  of  the.  manipulator. 

If  we  assume  f!  -  I  (identity  matrix)  and  N  =  0,  the  corresponding 
constant  coefficient  control  with  the  model  of  reference  is  written  according  to 
(5.2): 


u  —  —  L  f  - 
a  p 


i  E  +  U  . 

v  r 


(6.1) 


Using  the  equations  of  the  mathematical  model  of  the  manipulator  and  the 
reference  model  (4.3),  we  obtain  the  following  equation  for  the  error  e. 


e  =  M-1  u  -u  -  M-1N, 
a  r 

and  using  (6.1)  and  (6.2) 

e  -  -M‘1L  e-M^L  e+(M~1-I)u  -M_1N. 
p  v  r 


(6.2) 


(6.3) 


The  choice  of  large  gains  and  makes  the  term  -M  *LpC-M 


dominant  over  the  tern  (M  ^-Du^-M  as  soon  as  the  position  error  e  or  the 

velocity  error  e  becomes  large.  Then  the  equation  (6.3)  approximately 
reduces  to: 

c  =  -M  c-M  € .  (6.4) 

P  V 

Suppose  for  simplification  that  M  is  constant  or  that  its  relative  variation 
is  small  compared  with  that  of  e  (an  hypothesis  that  we  can  verify  later).  M 

being  positive  definite,  M  '  is  also  positive  definite  and  consequently  there 
exists  an  orthogonal  matrix  P  such  that: 

T  -1  -1 

PM  P  =  M  :  MD  is  a  diagonal  matrix  (6.5a) 

PTP  =  I.  (6.5b) 

If  we  let  e  =  PC  (6.6) 

and  make  use  of  (6.4)  and  6.5b)  we  obtain: 

**  T  -1  T  -1 

c  -  -  PM  L  P?  -  PM  L  PC  (6.7) 

p  v 

where  L  and  L  are  chosen  so  that: 

P  v 

L  =ir,L  =11  (1  >0,1  >0).  (6.8) 

p  p  v  v  p  V 

Lp  and  can  permute  with  the  matrix  P  in  the  equation  (6.7),  and 
using  (6.5a)  we  obtain: 

C=-L  C-LC  (6.9) 

p  v 

with  L'  =1  M D_1  ,  L'  =  1  M  :  positive  diagonal  matrices. 

The  relation  (6.9)  is  a  system  of  n  decoupled  equations  of  second  order, 
exponentially  stable.  Consequently  C  converges  to  0  exponentially  in  the 
same  way  as  e,  following  (6.6),  the  rate  of  convergence  being  the  more  rapid 
the  larger  the  gains  1^  and  1^. 

The  preceding  heuristic  reasoning  although  not  a  valid  demonstration, 
helps  us  justify  that  the  error  of  approximation  e  cannot  become  large.  If  an 
error  appears  at  any  instant,  then  it  diminishes  exponentially  and  more 
rapidly  as  the  gains  1^  and  1^  became  larger. 

We  now  examine  the  static  error  caused  by  the  gravity  termjn  N. 
Suppose  that  x  and  xr  converge  toward  constant  values.  Then  x  =  0, 

x  =  0,  u  =  0  (in  view  of  (4.3) J,  N(x,*x)  is  reduced  to  -u(x)  [in  view  of 
r  r  g 

(2.4))  if  we  suppose  the  friction  term  is  of  the  form  u  (x)  =  -k^x),  and 

equation  (6.3)  becomes  9 
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0 


(6.4) 


=  -  M-1  L  e  .  .  +  M“1u  (x) , 

p  stat  g 


from  which 


u  (x) 

c  =  — 2 — " - 

stat  . 


(6.5) 


showing  that  the  static  error  due  to  the  gravity  term  u  (x)  is  inversely 
proportional  to  the  gain  1^.  " 

Numerous  simulations  given  by  Ibarra-Zannatha  confirm  the  robustness  of 
large  gain  controls.  The  robustness  property  will  be  studied  further. 


Discretization  of  the  Control  Problem 

Using  relation  (2.1)  the  equation  of  the  mathematical  model  of  the  manip¬ 
ulator  can  also  be  written  in  the  form  of  the  equation  of  state: 

X  =  AX  +  B  (x)u  +  C (x,x)  (7.1) 

a 

with 


Let  us  assume  ug  is  numerically  computed  and  maintained  constant  during 

a  sampling  step  of  length  h.  In  integrating  (7.1)  over  a  sampling  step  we 
obtain: 


with 


X  (t+h) =AdX (t) +Bdua (t) +Cd 


Ad  =  e 


(7.2) 


I 

hi 

n 

n 

0 

I 

n 

B  • 


/t+h 
t 


eA(t+h-s)  B(x(s))ds  ;  and 


=  yt+heA(t+h-s)  c(x(s)>  £(s)]  ds> 
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r 


In  general,  it  is  not  possible  to  calculate  exactly  B^  and  C^. 


However, 


if  the  sampling  step  h  is  chosen  sufficiently  small,  we  can  assume  that  x  and 
x  vary  little  over  this  sampling  step,  as  do  B(x)  and  C(x,x).  B^  and 
can  then  be  approximated  as  follows: 


*t+h  A  (t+h-s) 


Bfx(t)] 


-(!•;>" 


(x(t)]  and 


(7.3a) 


t+h  A (t+h-s) 
e  dt 

t 


C  [x  (t)  ,x  (t)  ]  = 


'  (£)  M'1 


[X  (t) ] N [x  (t)  ,X(t)  ]  . 


(7.3b) 


From  equations  (7.2)  and  (7.3),  we  find  that  the  discretization  of  the 
problem  of  control  of  a  manipulator  leads  to  replacing  the  system  of  contin¬ 
uous  equations  (2.5)  by  the  discrete  system: 


X  (t+h)  =  AdX(t)  + 


uai  (t)  =  M  [x  (t)  ]  u'  (t) 
Ua2  =  N  (x  (t)  ,  x  (t)  ] 

ua(t)  =  ual(t)  +  ua2(t) 


(7.4a) 

(7.4b) 
(7.4c) 
(7. 4d) 


The  similarity  between  this  system  and  the  continuous  system  (2.5) 
shows  that  the  approach  we  have  considered  in  the  continuous  case  also 
applies  in  the  discrete  case.  The  system  shows  in  effect  that  the  problem  of 
computing  the  control  ug  is  decomposed  into: 

1.  A  control  problem  of  a  linear  decoupled  and  discretized  system 
(7.4a) . 

2.  A  modelling  problem  connected  with  the  relations  (7.4b)  and  (7.4c). 

In  particular,  the  property  of  robustness  of  control  with  large  gains  is 
preserved.  We  note  also  that  the  previous  discretization  eliminates  the  prob¬ 
lem  of  determining  the  second  derivative  of  x. 

Computation  of  the  Control  in  the  Frame  of  Reference  of  the  Operational 


Coordinates 

In  many  applications,  the  user  is  interested  chiefly  in  the  control  of  the 
terminal  arm  of  the  manipulator.  The  position  and  speed  of  the  arm  are 
usually  characterized  by  coordinates  different  from  the  manipulator's  gener¬ 
alized  coordinates,  which  are  less  practical  to  use.  General  coordinates 


10 


include,  for  example,  the  three  coordinates  in  a  fixed  frame  of  reference  of 
the  extremity  of  the  terminal  arm,  with  the  three  Euler  angles  in  this  frame 
of  reference  characterizing  the  inclination  of  the  terminal  arm.  To  solve  the 
control  problem  we  have  two  possibilities: 

1 .  Recompute  the  reference  trajectory  in  the  frame  of  reference  of 
generalized  coordinates,  then  compute  the  control  u  as  before.  The  biggest 

3 

difficulty  with  this  method  is  in  passing  from  the  operational  coordinates  to 
the  generalized  coordinates,  the  inverse  operation  being  generally  much 
easier. 

2.  Compute  directly  the  control  ug  from  the  operational  coordinates. 

The  possibility  returns  to  the  recomputation  of  the  model  of  the  manipulator 
in  the  frame  of  reference  of  operational  coordinates. 

We  show  that  for  the  second  possibility,  the  approach  developed  earlier 
in  this  report  remains  applicable. 

Let  y  be  the  vector  of  operational  coordinates  corresponding  to  the 
vector  x  of  the  generalized  coordinates  of  the  manipulator,  and  let  f(.)  be 
the  function  assumed  with  continuous  derivative  everywhere,  which  allows 
passing  from  x  to  y. 


y  =  f  (x) 

Let  J(x)  be  the  Jacobian  matrix  of  f  defined  as 

*«*>  -  {£• 

Following  (8.1)  and  (8.2) 
y  =  J  (x)  x. 


(8.1) 


(8.2) 


(8.3) 


In  taking  the  derivative  of  (8.3)  we  obtain: 
y  =  J  (x)  x  +  W  (x,x) 


(8.4) 


with 


W  (X/X) 


*T  1 

X  W  (x)  x 


x^V1  (x)  X 


(8.5) 


where  W  is  the  derivative  of  the  i  line  of  J  with  respect  to  x.  We  suppose 
for  simplification  that  the  matrix  J(x)  has  an  inverse.  (If  J  has  no  inver-e 
we  can  use  a  pseudo-inverse  of  the  matrix.) 

Following  (8.4)  we  have 


=  J_1(x) 


"  -1 
y-  J  (x) 


W  (x,x) 


(8.5) 
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Recalling  equation  (2.3'  of  the  mathematical  model  of  the  manipulator 


M  (x)  x  +  N(x,x)  =  u  .  (8.6) 

cl 

Combining  (8.5)  and  (8.6)  we  obtain 

M’ (x)  y  +  N ' (x, x)  =  (8.7a) 

u  =  JT(x)  (8.7b) 

a.  a 

with 

M' (x)  =  J-T(x)  M(x)  J-1(x)  and  (8.7c) 

N'  (x,x)=-M'  (x)W(x,x)  +  J  T(x)  N(x,x).  (8.7d) 


Comparison  of  (8.7a)  and  (8.6)  shows  that  the  two  relations  have  the 

same  form.  For  instance,  the  matrix  of  inertia  M(x)  is  replaced  by  the 

matrix  M'(x),  which  is  also  positive  definite  if  the  Jacobian  J(x)  has  an 

inverse.  It  follows  that  the  techniques  of  computing  the  control  u  (shown 

3 

previously)  are  transposable  to  the  computation  of  the  control  V  using  the 

operational  coordinates  of  the  manipulator.  The  control  u^  that  interests  us 

is  then  computed  from  according  to  relation  (8.7b).  The  idea  is  found  in 

M.  Takegaki  and  S.  Arimoto,  "An  Adaptive  Trajectory  Control  of  Manipula¬ 
tors,"  International  Journal  of  Control,  Vol  34,  No.  2  (1981),  pp  219-230. 

Conclusion 

This  report  has  presented  a  synthetic  approach  for  calculating  the 
control  of  manipulators.  The  initial  control  problem  is  broken  down  into 
linear  control  and  modelling  problems.  The  approach  allows  derivation  of 
numerous  schemes  (adaptive  or  not)  of  control  proposed  in  the  literature  and 
suggests  new  schemes.  It  has  been  shown  that  the  problem  of  modelling  is 
difficult  but  is  less  crucial  if  one  can  synthesize  robust  controls  that  are  not 
sensitive  to  errors  of  modelling.  A  heuristic  study,  which  will  later  be 
refined,  shows  that  large  gain  controls  are  robust  and  that  controls  with  a 
reference  model  are  naturally  conditioned  to  function  with  large  gains. 

To  give  the  study  its  correct  weight,  one  must  take  a  larger  view  and 
consider  the  control  problem  in  a  more  general  context  than  we  have  so  far. 
Some  difficulties  do  not  appear  with  our  mathematical  model  of  the  manipu¬ 
lator.  The  model  is  only  a  simplistic,  (although  very  useful)  representation  of 
the  physical  system  of  a  manipulator.  The  model  does  not  take  into  account, 
for  example,  the  flexibility  of  the  manipulator  arms  and  the  transmission 
organs.  Vile  have  not  considered  possible  errors  in  measuring  x  and  x,  which 
could  distort  in  an  important  way  the  performance  of  the  control.  It  is  easily 
conceivable  that  the  set  of  the  nonmodelled  terms,  which  are  fundamentally 


different  from  those  we  have  included,  could  change  the  basic  structure  of 
the  system  (in  the  case  of  the  flexibilities  when  they  are  important),  and 
make  risky  the  control  (in  the  case  of  noisy  measurements).  Such  terms  are 
not  modelled  now  in  studies  of  adaptive  control  for  manipulators,  but  they  can 
make  the  adaptation  of  control  gains  desirable  and  justify  the  study  of  more 
elaborate  controls. 

The  above  is,  in  our  view,  an  important  outcome  of  the  present  study, 
which  must  be  viewed  as  a  simple  element  of  analysis,  among  others,  of  the 
general  problem  of  the  control  of  manipulators. 
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